Atom Fock state preparation by trap reduction 
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We study the production of low atom number Fock states by reducing suddenly the potential trap 
in a ID strongly interacting (Tonks-Girardeau) gas. The fidelity of the Fock state preparation is 
characterized by the average and variance of the number of trapped atoms. Two different methods 
are considered: making the trap shallower (atom culling [A. M. Dudarev et al, Phys. Rev. Lett. 98, 
063001 (2007)], also termed "trap weakening" here) and making the trap narrower (trap squeezing). 
When used independently, the efficiency of both procedures is limited as a result of the truncation 
of the final state in momentum or position space with respect to the ideal atom number state. 
However, their combination provides a robust and efficient strategy to create ideal Fock states. 

PACS numbers: 32.80.Pj, 05.30.Jp, 03.75.Kk 



The generation of Fock states with a definite, con- 
trolled atomic number is a highly desirable objective both 
from fundamental and applied points of view. They may 
be useful for studying few particle interacting systems 
P, 0], entanglement [3], or number and spin-squeezed 
atomic systems d, 0] . Production of fotonic Fock states 
and interferometric schemes for sub-shot-noise sensi- 
tivity approaching the Heisenberg limit Q do also require 
input Fock states. 

A necessary step towards this goal is the development 
of atom counting devices paving the way into quantum 
atom statistics S]. Indeed a technique with nearly unit 
efficiency has already been demonstrated 19!]. Moreover, 
the recently proposed method of atom culling by mak- 
ing the atom trap shallower (in the following also termed 
"trap weakening") has achieved sub-Poissonian 

atom-number fiuctuations for 60 — 300 trapped atoms 
when adiabatic conditions were fulfilled, i.e., when the 
trap depth is varied slowly. (The reference case of Pois- 
sonian statistics is realized by the number of particles in 
a small volume of a classical ideal gas.) In this method 
the initial state is assumed to be a ground state for an 
unknown number of bosons, in general smaller than the 
maximum capacity of the initial trap (this is, the max- 
imum number of particles that can be confined in the 
trap). This capacity depends on the trap characteristics 
and interatomic interaction. As the barrier height of the 
trap is slowly reduced and the maximum capacity is sur- 
passed, the excess of atoms will leave the trap to produce 
eventually the Fock state corresponding to the maximum 
capacity of the final trap configuration. For a pictorial 
representation see Fig. 1 (upper panel). 

A theoretical analysis has shown the basic properties 
of atom culling regarding the final average number as a 
function of the trap well depth and atom-atom interac- 
tion strength, covering the limits of the Tonks-Girardeau 
(TG) gas and the mean-field regime a weaker de- 
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pendence on laser fluctuations of the height of the bar- 
rier that forms the trap will be favoured by strong in- 
teratomic interactions that separate the energy levels. 
This motivates the present work, in which we focus on 
the strongly interacting ID TG limit, optimal for atom 
culling, and examine the average number and fiuctua- 
tions of the trapped atoms. Particular attention is paid 
to the sudden regime, corresponding to the "worst case 
scenario" of an abrupt change from the initial to the final 
trap. We show that, even in this case, a state arbitrarily 
close to the ideal Fock state may be robustly produced 
by combining weakening and squeezing of the trap, the 
two basic processes represented schematically in Fig. 1. 

I. ATOM STATISTICS 

Ultracold bosonic atoms in waveguides tight enough 
so that the transverse degrees of freedom are frozen out, 
are well described by the Lieb-Liniger (LL) model [HI. 
In the strongly interacting limit jl2l . [l3| (for low densi- 
ties and/or large one-dimensional scattering length) a LL 
gas tends to the TG gas, which plays a distinguished role 
in atom statistics since its spatial antibunching has been 
predicted and observed [l^. The system exhibits 
"fermionization" [l^, i.e., the TG gas and its "dual" 
system of spin-polarized ideal fermions behave similarly, 
and share the same one-particle spatial density as well as 
any other local-correlation function, while differ on the 
non-local correlations. 

The fermionic many-body ground state wavefunc- 
tion of the dual system is built at time t = as a 
Slater determinant for Ni particles, iPf{xi, . . . ,X]\[.) = 
■^^=fdet^^f.^-^^ipl^{xk), where (Pn{x) is the n-th eigen- 
state of the initial trap, whose time evolution will 
be denoted by ipn{x,t) whenever the external trap 
is modified. The bosonic wave function, sym- 
metric under permutation of particles, is obtained 
from ipF by the Fermi-Bose (FB) mapping [l^ [l^ 
V'(xi, . . . ,XArJ = A{xi,. . . ,XNi)'^I^F{xi,. ■ ■ ,XN,), where 
= Y\i<j<k<N,^S^(^k - Xj) is the "antisymmet- 
ric unit function" . Since A does not include time 
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FIG. 1: (color online) Schematic potential change for trap 
weakening and squeezing, where Li and Lf are the initial 
and final well widths, and Vi and V/ the initial and final 
trap depths, respectively. Each solid line in the contour plot 
represents an isospectral family of traps at the threshold of 
a new bound state (for the initial trap 2mViLl/h^ = 25iv^, 
Cf being the final trap capacity). A transition between two 
given families is achieved reducing the trap capacity. The 
bound energy levels are pushed up and the higher ones cross 
the continuum threshold so that the excess of atoms escapes 
from the trap. 



reduced trap we need to calculate the atom-number fluc- 
tuations. We proceed by characterizing the TG trapped 
state by means of its variance tr^ = {N^{t)) — {N{t))^. 
First note the general relation 

{n{x)n{x'))t = (: n{x)n{x') :)< + d{x - x'){n{x))t, (1) 

where the number field operator is n{x) =^ {x)'^' (x) , 
^(x), ^'^(a;) are the annihilation and creation operators 
at point X, and : : denotes normal ordering. In partic- 
ular, within the Tonks- Girardeau regime. 



{n{x)n{x'))t = D{x,x';t) + 6{x - x')p{x,t), (2) 



with \2C 



D{x,x'-t) ^ N,{N, - 1) l[[dx,\^{x,x' ,X3, 



p{x,t)p{x',t) - |AA,,(a;,a;',t)|^ 



and 



Ni 



ANi{x,x';t) = ^ (pn{x,t)*(pn{x',t). 



(3) 
(4) 



n=l 



The mean value of the number of particles within the 
trap and of its square can be obtained by integrating 



{N{t)) = 



dxp{x, t), 



dxdx' {n{x)n{x'))t, (5) 



where i + ^ is large enough to include the bound- 

state tails in coordinate space. (For the trap configura- 
tion of Fig. 1 each bound state has a penetration length 
{^j} beyond the well width i, so ^ = maxj^j.) The atom 
number variance reads finally 



a%{t) = (Nit)) 



dxdx' An.{x, x' ,t). (6) 



From it one can infer Poissonian statistics if a'j^/ {N) > 1 
and sub-Poissonian as long as ajf/{N) < 1. For a Fock 
state (Tat = 0. 



explicitly, it is also valid when the trap Hamilto- 
nian is altered, and the time-dependent density pro- 
file resulting from this change can be calculated as 
[11] p{x,t) = NiJ \^l;{x,X2, ■ ■ ■ ,XN,;t)\'^dx2 ■ ■ ■ dxN^ = 
J2n=i\'Pn{x,t)\'^ . By reducing the trap capacity some 
of the Ni atoms initially confined may escape and N 
will remain trapped. To determine whether or not sub- 
Poissonian statistics or a Fock state are achieved in the 



II. THE SUDDEN APPROXIMATION 

The requirement of adiabaticity, i.e., of a slow trap 
change, is a handicap which one would like to overcome. 
Achieving good fidelity with respect to the desired Fock 
state may require exceedingly long times, a fact that is 
even more critical whenever the interactions are finite, 
this is, for the Lieb-Liniger gas in which the splitting be- 
tween adjacent levels (Bethe roots) diminishes. It is thus 
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useful to examine the opposite limit corresponding to a 
sudden trap change. For the TG gas, we shall find gen- 
eral and exact results which are a useful guide since small 
deviations from the sudden limit increasing the switching 
time may only improve the fidelity. 



In what follows we shall thus discuss the preparation 
of Fock states by an abrupt change of the trap poten- 
tial to reduce its capacity. Even though the arguments 
and results of this section are rather general, consider for 
concreteness the simple square trap configuration of Fig. 
1, with an infinite wall on one side and a flat potential 
(zero potential energy) on the other side, 

{oo , X < 
, 0<x<L^ , (7) 
, X > La 

where the subindex a — i, f refers to the initial and final 
configuration respectively. The corresponding eigenvalue 
problem is solved in the appendix, where both bound 
and scattering states are described. At t = the trap 
with shape Vi{x), which holds an unknown number of 
particles Ni (lower or equal than the initial capacity Ci) is 
suddenly modified to the final trap Vf (x) , which supports 
Cf bound states ifij{x), j = 1, ■■■,Cf, of energy < 0. 
The process may consist on "weakening" the trap {Vf < 
Vi, Lf = Lj), "squeezing" it {Lf < Li, Vf = V), or a 
combination of the two {Vf < Vi and Lf < Li). 



The continuum part of the spectrum of the one-particle 
Hamiltonian with potential Vf {x) is spanned by the scat- 
tering states x{(a;), labeled by the incident wavenumber 
k. It follows from standard scattering theory [l^, using 
the Riemann-Lebesgue lemma ^2l\ , that the contribution 
of the continuum states in the trap region {x < L+) van- 
ishes asymptotically as t oo, 

+ / dk{x\xi){xiW)^)e'''^ 
Jo 

Cf 

~ J2{^\^i){vi\'Pn)e-'''''/\ (8) 

so that the dynamics in the trap is finally governed by the 
discrete part of the spectrum {ipj\j = 1, . . . , C/}. There- 
fore, asymptotically, the mean number and variance of 
trapped atoms are 

Cf 

(iV(oo)) = ^(^:jA;|^j,)=X^(^^/|A,|^J), (9) 
ji=i j=i 



note that {N{oo)) < Cf, and 

a^(oo) = (^(oo))- ^ MJAf\^l,)\' 

n,m— 1 

Cf 

= ^(^J|(A,-A,A/A,)|^J), (10) 

3 

where 

Cf 

A/ = El'^J)(^Jl (11) 

is the projector onto the final bound states and 

A,^Y.\^n){Vu\ (12) 

ri=l 

the projector onto the bound states occupied by the ini- 
tial state. We may thus conclude that trap reduction 
can actually lead to the creation of Fock states with 
(iV(oo)) = Cf and ajf{oo) = when the initial states 
span the final ones, 

A/CA,. (13) 

A time scale for the validity of the asymptotic regime, 
after the trap switch, is provided by the lifetime of the 
lowest resonance of the final trap. A simple semiclassical 
estimate is r = Lf{m/2Vfy^^, assuming the escape of 
a classical particle from the well and approximating the 
resonance kinetic energy by the potential depth. To pre- 
pare a ^^Na Fock state in final trap of with L/ ~ 50 fim, 
supporting Nf = 10 bound states, the asymptotic regime 
is approached for i > r ~ 30 ms. We insist though, that 
any slower potential change will play in favour of the fi- 
delity of the resulting Fock state until the time scale in 
which losses and decoherence begin to play a role. 



III. TRAP WEAKENING 

A good guidance for Fock state preparation is provided 
by Eq. PH)) . which leads to the requirement for A^ to be 
an extension of A/. This result is model- independent, 
and in particular it holds irrespective of the smoothness 
of the trapping potential. For illustration purposes we 
shall consider the square trap in Eq. ^ shown in Fig. [T] 

In a trap- weakening scheme the potential is made shal- 
lower by reducing the depth from an initial value Vi to 
Vf (while Li — Lf, see Fig. [IJ, a procedure which 
has been successfully implemented to prepare states with 
sub-Poissonian statistics [l3|- Though, in practice, the 
value of Ni (and Cf) cannot be arbitrarily large be- 
cause the Tonks-Girardeau regime requires a linear den- 



4 



0.94 r 
,0.90^. 

0.86 r 
0.82^ 



weakening 



a) : 



200 



0.92 



S 



400 600 

I ' I 

squeezing 



800 



b) 



^ 0.91^ 
0.90- 



20 



40 



c. 



60 



80 



FIG. 2: (color online) Limited efficiency on pure squeezing 
and weakening techniques. For a completely loaded initial 
trap {Ni = Ci) of depth Ui = Cfn^, the resolution of a given 
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{n = 5, Uf = 257r^), 



bound state of the final trap, 
quantified by the measure F„ — (if(^\Ai\ipl^), and shown to be 
limited in a) pure squeezing, and b) pure weakening. For in- 
creasing capacity of the initial trap, d, F„ tends to the prob- 
ability to find the particle in the well region [0, L^] (weakening 
case, solid line, see Eq. (|14|) ). or to the probability to find the 
particle in a momentum space region (squeezing case, solid 
line, see Ea. (fT5)l '). 
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FIG. 3: (color online) Asymptotic mean value (up) and vari- 
ance (down) of the atom number distribution of a Tonks- 
Girardeau gas obtained by sudden trap reduction as a func- 
tion of the width ratio between the final and initial trap. 
The initial trap with Ui/n'^ — 10* supports a maximum 
of Ci — 100 bound states, while the final configuration 
Uf/iv'^ = 10^ is limited to Cf = 10. Different filling factors 
are considered for the initial trap. The left and right edges 
of each curve correspond to pure squeezing and weakening 
respectively, the later being remarkably less sensitive to the 
trap filling as shown in the inset. Any other point combines 
weakening and squeezing. 



sity n ~ 1 , it is useful to consider a large number 

of particles in an initial box- like trap, Ni oo, for which 
Ai becomes the projector in the interval [0, L;], 

-X[o,L.](^). (14) 

This asymptotic behavior is depicted in Fig. [2^. Prepa- 
ration of ideal atom number states by trap weakening is 
thus hindered by the suppression of the coordinate space 
tails which leak beyond the well along a given penetration 
length for each ip^ [25l |. 



IV. TRAP SQUEEZING 

There is a simple alternative to trap weakening to 
achieve high-fidelity Fock states: atom-trap squeezing. 
Starting with a state of an unknown number of particles 
Ni, the trap width is reduced from an initial value Li 
to Lf keeping the depth constant {Vi — Vf), as shown 
in Fig. [1] (middle panel). The final trap supports Cf 
bound states so that the excess of atoms is squeezed 
out of the trap. From a comparison of initial and fi- 
nal energy levels, it is clear that a minimum number of 
initial particles is required for trap squeezing to work. 
Using for an estimate the levels of the infinite well we 
get Ni > CfLi/Lf. Trap squeezing works optimally for 
initial traps filled with atoms to the brim but it is not 
robust against partial filling. It is also less sensitive to 
threshold effects than trap weakening, in particular for 
low atom numbers. 



For a wide, filled initial trap, 

Ai~ / dk\k+){k+\ (15) 
Jo 

as L — > oo, where {x\k~^) = ^2/7rsin kx, x > 0, satisfy- 
ing = S{k—k') and the cut-off is at k ~ U^^^/Li, 
in terms of the dimensionless parameter Ui = 2'mLjVi / h? 
{Uf is defined similarly in terms of the final values). This 
is illustrated in Fig. ^jp. Therefore, trap squeezing may 
limit the fidelity of the final Fock state preparation due 
to the truncation of the tails "in momentum space" , in 
the sense of Eq. HID, for fc > Ul^^/L,. 



V. COMBINED WEAKENING AND 
SQUEEZING 

From the previous discussion it follows that the opti- 
mal potential trap change to fuUfiU Af C A^ is a com- 
bination of weakening and squeezing. Let us choose two 
different families of isospectral traps characterized by Ui 
and Uf, supporting Ci and Cf bound states respectively 
(In general Ci > Ni because of partial filling of the ini- 
tial trap, the filling factor being the ratio Ni/Ci). The 
energy level of the highest occupied state measured from 
the bottom of the trap, analogous to the Fermi level in 
the dual system, is denoted by e^, which depends on the 
filling factor of the JJi-trap, while for the ideal final Fock 
state £/ ~ V/. 

From an arbitrary potential (Li,Vi) in the [/^-family. 
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FIG. 4: (color online) Asymptotic mean value (up) and vari- 
ance (down) of the atom number distribution of a Tonks- 
Girardeau gas obtained by a sudden change of the trap po- 
tential for Ui/iv'^ = 10" (solid line), 2.25 x lO" (dotted line) 
and 4 x 10^ (dashed line) keeping a constant initial filling fac- 
tor of 90%. Uf/ix^ = 10^ with Cf = 10. The high-efficiency 
region increases with Ui reducing the critical ratio between 
the final and initial trap widths L'^^/Li. 



the efficiency of the weakening-squeezing combination 
leading to a potential of the Uj family varies with the 
filling factor and with the ratio of widths in the final 
and initial trap Lf/Li. Figure [3] shows both the mean 
atom number and variance of the final states for differ- 
ent ratios and preparation states. Below a critical final 
width — {Uf/UiY/'^Li, the physical final depth Vf 
is larger than the initial one Vi and we disregard this 
possibility since the efficiency is very poor as expected 
from the failure of the condition (|13p for £f > e^. The 
ratio Lf/Li = 1 is the limit of pure trap weakening, and 
/ Li that of pure trap squeezing; the limited efficiency 
of both extreme cases can be noticed in Fig. [3l However, 
the combined process achieves pure Fock states and is 
robust with respect to different fillings in the initial trap 
and for a wide range of final configurations. Moreover, 
the range of final configurations for which high-fidelity 
states are obtained increases with the initial potential Ui 
keeping the filling constant as shown in Fig. [H Hence, 
a recipe to create a Fock state |C/) in a trap of width 
Lf and depth Vf will be as follows: choose the width of 
the initial trap broad enough in the sense Li Lf + 
(r > 1), where ^ is the penetration length of the state 



, in the final trap. Then, make sure that Si > 



such a way that the final state is contained in the ini- 
tial subspace both in momentum and coordinate space. 
A deep and broad initial trap, with respect to the final 
one, provides in summary a safe starting point to create 
a Fock state by sudden (or otherwise) trap reduction. 



VI. DISCUSSION AND CONCLUSION 

In this work we have studied and compared strategies 
for atom Fock state creation in Tonks- Girardeau regime. 



We have shown that Fock states can be prepared even un- 
der a sudden trap-potential change. Thanks to the anal- 
ysis of the atom number variance, we have determined 
that the key condition for Fock state creation is that the 
initial occupied bound states span the space of the fi- 
nal ones. This holds regardless of the trap shape details, 
and in particular does not depend on the potential-trap 
model. A combination of trap weakening and squeezing 
allows to resolve the ideal Fock state both in momen- 
tum and coordinate space. We close by noting that the 
Tonks-Girardeau regime is optimal with respect to the 
strength of interactions. In this regime, the three-body 
correlation function 33(7) tends to vanish and therefore 
the losses of atoms from the trap by inelastic collisions 
are negligible ^28!|. For gases with finite interactions in 
tight- waveguides, when the Lieb-Liniger model holds, the 
quasimomenta obtained as solution of the Bethe equa- 
tions [33, [31I, [S^l , are closer to each other for weaker 
interatomic interactions, whence it follows that the re- 
quired time scale for the dynamics to be adiabatic is even 
larger than in the Tonks-Girardeau regime. Large spac- 
ings in the Bethe roots also imply that less precision is 
required in the control of Vf . Given that the interatomic 
interactions can be tuned through the Feschbach reso- 
nance technique, one can optimize the atom Fock state 
preparation by putting the system within a strong inter- 
action regime in a first stage, followed by the controlled 
reduction of the potential trap (weakening and squeez- 
ing) , and finally turning off slowly the interaction. 
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APPENDIX A 

In this appendix we describe the spectrum of the 
Hamiltonian for the potential in Eq. ([7]) considered in 
the numerical examples. Both the initial and final traps 
have the same functional dependence. Here we consider 
the general case for a trap of width L and depth V (drop- 
ping the index a = i, f for compactness), whose spectrum 
can be easily determined using matching conditions in 
the wavefunction and its derivatives. The trap supports 
a finite set of C bound states (its capacity) 



ifj(x) =J\fj 



0<x<L 



sm{qjx) 
sin(gfjL)e~''j(''~-^) , x> L 



(Al) 
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with j = 1, . . . ,C and normalization constant 

-1/2 



L sm{2qjL) ^ sin (qjL) 



2k,- 



(A2) 



The eigenvalues are Ej = [hqjY /2m— F < where {qj} 
satisfy the trascendental equation Kj = —qj cot qjL, with 
Kj = {2mV/h^ — q]y^^ > 0. The scattering states have 
the form 



Xk{x) 



1 



C A sing's; 
2iT I e-''^^ - S{k)e'' 



0<x<L 
x> L 



,(A3) 



where q = y^k"^ + 2mV/h? and the coefRcient A and the 
scattering matrix S{k) are determined by imposing the 
usual matching conditions, 



A 
S{k) 



2ike- 



ikL 



= e 



q cos qL — ik sin qL ' 

_2ikL 9 1^ + si"^ 9-^ 
g cos qL — ik sin qL 



(A4) 
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